PSU/TH/139 
[icp-th/9312242 
December 1993 

INSTANTON-INDUCED PARTICLE PRODUCTION 
IN DEEP INELASTIC SCATTERING [] 

IAN BALITSKY f\ 

Physics Department, Penn State Univ., 104 Davey Lab. 
University Park, PA 16802, USA 

ABSTRACT 

I discuss the instanton-induced contributions to the coefficient functions in front of 
parton densities. They correspond to the spherically symmetric particle production 
at the level of 1CP 2 — 10~ 5 of the total cross section of deep inelastic scattering from 
the parton. 

1. Introduction 

The deep inelastic lepton-hadron scattering at large momentum transfers Q 2 
and not too small values of the Bjorken scaling variable x = Q 2 /2pq is studied in much 
detail and presents a classical example for the application of perturbative QCD. The 
factorization theorems allow one to separate the Q 2 dependence of the structure 
functions in coefficient functions Ci(x,Q 2 / p 2 ,a s (p 2 )) in front of parton (quark and 
gluon) distributions of leading twist Pi(x,p 2 ,a s (p 2 )) 

F 2 (x, Q 2 ) = EiCi(x, Q 2 /p 2 , a s (p 2 )) <g> P t {x, p 2 , q,^ 2 )), (1) 

where the summation goes over all species of partons, and fi is the scale separating 
"hard" and "soft" contributions to the cross section. At p 2 = Q 2 the coefficient 
functions can be calculated perturbatively and are expanded in power series in the 
strong coupling 

C(x,l,a s (Q 2 )) = Co(x) + ^^C 1 (x)+ ( C 2 (x) + . . . (2) 

whereas their evolution with p 2 is given by famous Dokshitzer-Gribov-Lipatov- 
Altarelli-Parisi equations. Going over to a low normalization point p 2 ~ lGeV, 
one obtains the structure functions expressed in terms of the parton distributions 
in the nucleon at this reference scale. The parton distributions absorb all the in- 
formation about the dynamics of large distances and are fundamental quantities 
extracted from the experiment. Provided the parton distributions are known, all 
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the dependence of the structure functions on the momentum transfer is calculable 
and is contained in the coefficient functions Cj. Corrections to this simple picture 
come within perturbation theory from the parton distributions of higher twists and 
are suppressed by powers of the large momentum Q 2 . 

The picture described above presents a part of the common wisdom about 
hard processes in the QCD, and in a more or less detailed presentation can be found 
in any textbook. Less widely known is the fact that from the theoretical point of 
view this picture is not complete. An indication that some contributions may be 
missing comes from the asymptotic nature of the perturbative series in (||). This 
series is non-Borel-summable, which means that any attempt to attribute a quanti- 
tative meaning to the sum of the series in (0) would produce an exponentially small 
imaginary part ~ iexp{-const • 7r/a s (Q 2 )}, which is to be cancelled by the imaginary 
part coming from nonperturbative contributions. Thus, separation between pertur- 
bative and nonperturbative pieces in the cross section as the ones which contribute 
to the coefficient function and the parton distribution, respectively, cannot be rigor- 
ous. A modern discussion of the asymptotical properties of the perturbation series 
in QCD can be found in [0,0. 

In addition to imaginary exponential corrections which must cancel identically 
against the corresponding nonperturbative contributions, the coefficient functions 
may acquire also real exponential corrections, which potentially produce observ- 
able effects. In this talk I shall report on recent results by V.Braun and myself 
H |J, indicating that these corrections are indeed present. We have found that 
the deep inelastic cross section indeed possesses exponential contributions of the 
form F(x)exp[-4:TrS(x)/a s (Q 2 )], where S(x), F(x) are certain functions of Bjorken x, 
which we are able to calculate in a certain kinematical domain. Since the exper- 
imental data are becoming more and more precise, it is of acute interest to find 
a boundary for a possible accuracy of the perturbative approach, which is set by 
nonperturbative effects. Our study has been fuelled by recent findings of an en- 
hancement of instanton-induced effects at high energies in a related problem of the 
violation of baryon number in the electroweak theory [HP]. In the case of QCD the 
instanton-induced effects could turn out to be significant at high energies, despite 
the fact that they correspond formally to contributions of a very high fractional 
twist exp(-4nS(x)/a s (Q 2 )) ~ (A 2 QCD /Q 2 ) bs ^ . 

Note that, in principle, the instanton contributions in QCD are infrared- 
unstable - in a typical situation integrations over the instanton size are strongly 
IR-divergent. However, this problem does not affect calculation of instanton contri- 
butions to the coefficient functions which are IR-protected, as we shall see below 
(the detailed discussion can be found in H). 

2. Instanton contribution to the structure function of a gluon 

Let us start from the instanton-induced contribution to the coefficient func- 
tion in front of gluon density. Following Zakharov |7] , we calculate the cross section 
of the 7*5 scattering by means of the optical theorem. The trick is to evaluate the 



Figure 1: The instanton- induced contribution to the structure function of a gluon (b) and of a quark (c,d). 
Wavy lines are (nonperturbative) gluons. Solid lines are quark zero modes in the case that they are ending 
at the instanton (antiinstanton), and quark propagators at the II background otherwise. 



contribution to the functional integral coming from the vicinity of the instanton- 
antiinstanton configuration in Euclidian space, and calculate the cross section by 
the analytical continuation to Minkowski space and by taking the imaginary part. 
The true small parameter in our calculation is the value of the coupling constant 
at the scale Q 2 , which ensures that the effective instanton size is sufficiently small. 
We prefer to start, however, from the well-studied situation where not only the 
instanton sizes are small, but also the II separation is much larger than these sizes. 
As we shall see below, this corresponds to 1 - x < 1 (but of course 1 - x > a s still) . 
After that we shall move to smaller x (which correspond to strongly interacting / 
and I) trying to be accurate to collect to the semiclassical accuracy all the depen- 
dence on p 2 /R 2 in the exponent. To this end, we shall have in mind the valley method 



IT] , in which all the dependence on the II separation is absorbed in the action 
on the II configuration. However, we do not take into account corrections of 
order p 2 /R 2 in the preexponent, and to this accuracy need the first nontrivial term 
only in the cluster expansion of the quark propagator at the II background ||: 
(x|V7/V jr |0) =/dz(x|Vr 2 V 1 |z)^ s |(z|V 2 V- 2 |0). 

The leading contribution to the gluon matrix element of the T-product of 
the electromagnetic currents is given by the following expression 

dze«*{A a (p), \\T{^{z) ]v {0)}\A a {p), X} 11 = 

= Y, e lJ du J / y- d fa) J dR J dT J dzeiqz 



X 
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~ liDV^o p 4 e^Tr e -S 
x ( a t a )»/~i { a o (O)a l ,(O|V 2 - 2 V 2 aViVr 2 |z)^0o^) 

+a t K O (z)CT M (z|Vr 2 VidV 2 V 2 - 2 |0)CT I ,Ko(0) + (fi <-> v, z <-> 0) + . . .} (3) 

which corresponds to the diagram shown in Fig. la. The full expression contains 
many more terms H which are not shown because we have found that all of them 



are of order 0(a s (Q 2 )) compared to the expression in Eq.|3|. The subscript '1' refers 
to the antiinstanton with the size pi and the position of the center xj = R + T, and 
the subscript '2' refers to the instanton with the size p 2 and the center at xj = T. 
We use conventional notations V = V M o- M and V = V M d- M , etc, where a™ 61 = (—ia, 1), 
o>d Q = (+i(T, 1), and a are the standard Pauli matrices. Also, we write the quark zero 

modes in terms of the two-component Weil spinors tpo = y ^° J , V4 = (<t>o «o), and a 

and a* denote the overlap integrals a = - j dx (Rdn), = - J dx {4>B(f>). Here Sjj is the 
action of the instanton-antiinstanton configuration and £ = R +Pl+P2 is the conformal 

° S PlP2 

parameter []10j (the normalization is such as S(£) = 1 for an infinitely separated 
instanton and antiinstanton). Writing the action as a function of £ ensures that 
the interaction between instantons is small in two different limits: for a widely 
separated II pair, and for a small instanton put inside a big (anti) instanton, which 
are related to each other by the conformal transformation. In the limit of large £ 
the expansion of for the dominating maximum attractive II orientation reads 
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s(0 = (i-^ + o(H0/^)) (4) 

where the l/£ 2 term corresponds to a slightly corrected dipole-dipole interaction. 
Thus, the action S(£) decreases with the distance between instantons, so that the 
instanton and the instanton effectively attract each other. This attraction results 
in the exponential increase of the cross section — the effect found by Ringwald || . 
Further terms in the expansion of the action can be obtained by the so-called val- 
ley method [11|, and a typical solution (conformal valley) |I(J gives a monotonous 



function of the conformal parameter, which turns to zero at R — ► 0. In the tradi- 
tional language, the valley approach corresponds to the summation of all so called 
"soft-soft" corrections arising from the particle interaction in the final state. Main 
problem is in the evaluation of "hard- hard" corrections [|12|, which come from parti- 
cle interaction in the initial state. These corrections are likely to decrease the cross 
section, and in physical terms must take into account an (exponentially small) 
overlap between the initial state, which involves a few hard quanta, with the semi- 



classical final state [[13]]. Thus, the instanton-antiinstanton action is substituted by 
an effective "holy grail" function, which determines the leading exponential factor 
for the semiclassical production at high energies, and which received a lot of atten- 



tion in recent years. Unitarity arguments JTJ, |15| suggest that the decrease of the 
action will stop at values of order S(£) ~ 0.5. In a recent preprint [|16| Diakonov and 
Petrov argue that S(g) indeed decreases up to the value 1/2 at a certain energy of 
order of the sphaleron mass, and then starts to increase, so that the semiclassical 
production cross section is resonance-like. The question seems to us to be not set- 
tled finally. In this study, we have taken the value S = 1/2 as a reasonable guess 
for the residual suppression, and assumed that the behavior of the "true" function 
S(g) for S(£) > 1/2 is close to that given by the conformal valley [IIJ. The latter 
assumption is supported by numerical studies, e.g. in [T^]. 



We shall see below that the leading contribution in the strong coupling comes 



from the following regions of integration: 

z 2 ~ l/(Q 2 a s ) 
{z-R-Tf+pl ~ T 2 + p\ ~ z 2 /a s 

{z-R-T) 2 ~ T 2 - R 2 ~ p 2 - p 2 - z 2 /a s (5) 

and additionally p 2 /R 2 ~ 1 — a; when x is close to 1. Note that these regions of 
integration correspond to imaginary part of the II contribution so effectively the z? 
are negative. 

Since z 2 is small we can use the lightcone expansion (see e.g. for the 

quark propagator in the II background. Using the explicit expressions for the prop- 



agators from ]20] we find: 

Ro(x)a^(x\X7- 2 V 1 dV 2 X7- 2 \0)<j„K Q {0) = ( 
1 f\ {P1P2? 12 1 



2^ J Q ai [{ Z -R-T) 2 + p 2 ] 2 [T 2 + p 2 ] 2 ^W^\ Z 4 W 



, (jz-R-T) 
Pl (jz-R-Tf 
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' y/1 + pf/ijz - R-T) 2 d 1 + p y( lz - T f 



2 {iz-T) 
P2 (l*-T) 2 



Omitted terms have turned out to be of order 0(a s ). 

Let us at first consider the toy example of the Eq.^| without extra integration 
over 7 which brings only technical complexities: 

d p\, f d p\(„2^ [ jd *p*-£a-A) r(mi)r(n) r(m 2 ) 



^ rfT ^(P?)" 1 ^(pS)"" dRe ~™s 



Z 2n J J p\ 1 7 Pl V [(z-J?-r) 2 +p?] m i [T 2 +p 2 ] m 2 

(7) 

This integral diverges at p — ► oo. However, the divergent part corresponding 
to instanton with size p — Aq C d contributes only to parton densities and not to the 
coefficients in front of them. Moreover, these divergent parts possess no imaginary 
part so we shall imply them subtracted in what follows. (Strictly speaking we need 
fi-m subtractions of the type ((x — R — T) 2 + p 2 Y 2 - p~ 4 + 2(x — R — T) 2 p~ 6 + ...). After 
performing the integrations one obtains 

J {z 2 Y~ l J 2^+^-^(l+u)^+^- l[ a s (l + u) 2) {) 

where we use the notation I = //i + //2 — mi — mi + 4 and u = 1 — u. After continuation 
to Minkowski space the imaginary part of this integral take the form: 

R(n —I) r Hi+H2-2n+2l+2(i , ^n-l-fn-^ OA-n- 4,/, 6n 

fQ2\2-n+l 2 m i+ m 2+2n-2l~3 a s ^ 2 

where x = 1 - x i) = ^-fej^y and £ is now 2(1 + x)/x. It is easy to see now that 
the characteristic distances correspond to Eq|| 

One may also account for for the p dependence in the argument of a s . Careful 
analysis Q shows that ^ should be changed to + 2b ( where b = ^N c - |n/) and 
the argument of a s obeys the equation 

4tt 12(6,-2) 



Figure 2: The non-perturbative scale in deep inelastic scattering (instanton size p* 1 ), corresponding to 
the solution of equation Eqjl^ as a function of Q and for ~ 0.5 — 0.6 (f* ~ 3 — 4). 



A numerical solution of this equation for the particular expression of the action S(£) 
corresponding to the conformal instanton-antiinstanton valley is shown in Fig. 2. 
Note that the difference between the hard scale Q 2 and the effective scale for non- 
perturbative effects p~ 2 is numerically very large. This is a new situation compared 
to calculations of instanton-induced contributions to two-point correlation func- 

where the size of the instanton is of order of the large 
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tions, see e.g. 

virtuality. The effect is that the instanton-induced contributions to deep inelastic 
scattering may turn out to be non-negligible at the values Q 2 ~ lOOOGeV 2 , which are 
conventionally considered as a safe domain for perturbative QCD. 

Now, in order to find the 7*5 amplitude Eq.|3| we should take the real answer 
Eq.[7| instead of our toy example. After a considerable algebra (for details see || ) 
we obtain the following answer for the II contribution to the structure function of 
a real gluon: 
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>m 2 bs&)[bs&) - 1] 
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exp 
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(11) 



To our accuracy, we find that the instanton- induced contributions obey the Callan- 
Gross relation F 2 (x, Q 2 ) = 2xF 1 (x, Q 2 ). 

The expression in Eq.|Tl| presents our main result. It gives the exponential 
correction to the coefficient function in front of the gluon distribution of the lead- 
ing twist in Eq.|2|. The exponential factor is exact to the accuracy of Eq.[|. The 
preexponential factor is calculated to leading accuracy in the strong coupling and 
up to corrections of order 0(1 —x). The corresponding contribution to the structure 
function of the nucleon is obtained in a usual way, making a convolution of ( |TT| ) 
with a distribution of gluons in the proton at the scale pi. 



The instanton-antiinstanton contribution to the structure function of a quark 
contains a similar contribution shown in Fig. lb. The answer reads 



F$ q) (x, Q 2 ) 



E 
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128 rf 2 ^ 9 / 2 /16 

81i8 65(^)[^.)-l] U 3 

15/2 



exp 



47T 



26 5(6) 



(12) 



However, in this case additional contributions exist of the type shown in Fig.lc. 
They are finite (the integral over instanton size is cut off at p 2 — x 2 /a s ), but the 
relevant instanton-antiinstanton separation R is small, of order p. This probably 
means that the structure of nonperturbative contributions to quark distributions 
is more complicated. This question is under study. The answer given in Eq.O 



presents the contribution of the particular saddle point in Eq. 10 



3. Value of cross section and structure of final state for instanton-induced 
particle production 

The instanton-induced contribution to the structure function of a gluon in 



Eq.JT] is shown as a function of Bjorken x for different values of Q ~ 10 - lOOGeV 
in Fig. 3. The contribution of the box graph is shown by dots for comparison. 
The low boundary for possible values of Q is determined by the condition that 
the effective instanton size is not too large. At Q = lOGeV we find p* ~ lGeV" 1 , 
cf. Fig. 2. This value is sufficiently small, so that instantons are not distorted too 
strongly by large-scale vacuum fluctuations. Another limitation is that the valley 
approach to the calculation of the "holy grail" function S(£) is likely to be justified at 
S(0 > 1/2, which translates to the condition that x > 0.3-0.35. Numerical results are 
strongly sensitive to the particular value of the QCD scale parameter. We use the 
two-loop expression for the coupling with three active flavors, and the value Aj9 = 
365MeV which corresponds to the coupling at the scale of r-lepton mass a s (m T ) = 0.33 
2~I|j. Since the dependence on the coupling is exponential, the 20% increase of 



a s (p*) induces the increase of the cross section by almost an order of magnitude! 
Together with uncertainties in the function and in the preexponential factor, 
this indicates that the particular curves given in Fig. 3 should not be taken too 
seriously, and rather give a target for further theoretical (and experimental?) studies 
to shoot at. 

To summarize, we have found that instantons produce a well-defined and 
calculable contribution to the cross section of deep inelastic scattering for sufficiently 
large values of x and large Q 2 ~ 100 - lOOOGeF 2 , which turns out, however, to be 
rather small — of order 10~ 2 - 10~ 5 compared to the perturbative cross section. This 
means that the accuracy of standard perturbative analysis is sufficiently high, and 
that there is not much hope to observe the instanton-induced contributions to the 
total deep inelastic cross section experimentally. However, instantons are likely to 
produce events with a very specific structure of the final state, and such peculiarities 
may be subject to experimental search. The dominating Feynman diagrams in our 



Figure 3: Nonperturbative contribution to the structure function Fi(x,Q 2 ) of a real gluon ([H]) as a 
function of x for different values of Q (solid curves). The leading perturbative contribution is shown for 
comparison by dots. The dashed curves show lines with the constant effective value of the action on the II 
configuration. 

calculation correspond to 2n/a(p*) ~ 15 gluons and 2n/ -1 = 5 quarks in the final 
state with the low energy of order p^ 1 ~ 1 GeV. They are produced in the spherically 
symmetric way in the cm. frame of the partons colliding through the instanton.( 
The transverse momentum of the quark coming to the instanton is k\ ~ Q 2 a s ~ few 
Gev and so is the transverse momentum of the current quark jet) . It is not likely 
that quarks and gluons emitted from the instanton can be resolved as minijets (they 
have k± ~ Qa s ~ 1 Gev), and we rather expect a spherically symmetric production 
of final state hadrons in this frame. The effect is likely to be resonance-like, that is 
present in a narrow interval of values of Bjorken x of order 0.25-0.35 (in the parton- 
parton collision). In any case, finding of an instanton-induced particle production 
at high energies is a challenging problem, and further theoretical efforts are needed 
to put it as a practical proposal to experimentalists. 

5. Acknowledgements 

I would like to thank V.M.Braun for a very rewarding collaboration. It is a 



pleasure to acknowledge also useful discussions with J. C. Collins, E.M. Levin, A. H.Mueller, 
and M.I. Strikman.This work was supported by the US Department of Energy under 
the grant DE-FG02-90ER-40577. 

6. References 

1. A.H. Mueller, The QCD Perturbation Series, Talk given at workshop on QCD: 
20 Years Later, Aachen, 9-13 June 1992, CU-TP-573. 

2. V.I. Zakharov, Nucl. Phys. B385 (1992) 452. 

3. I.I. Balitsky and V.M. Braun, Phys. Rev. D47 (1993) 1879. 

4. I.I. Balitsky and V.M. Braun, Phys. Lett. B314 (1993) 237. 

5. A. Ringwald, Nucl. Phys. B330 (1990) 1. 

6. M.P. Mattis, Phys. Rep. 214 (1992) 159; 
P.G. Tinyakov, Preprint CERN-TH.6708/92. 

7. V.I. Zakharov, Preprint TPI-MINN-90/7-T (1990); Nucl. Phys. B371 (1992) 
637. 

8. N. Andrei and D.J. Gross, Phys. Rev. D18 (1978) 468. 

9. V.V. Khoze and A. Ringwald, Phys. Lett. B259 (1991) 106. 

10. A.V. Yung, Nucl. Phys. B297 (1988) 47. 

11. I.I. Balitsky and A.V. Yung, Phys. Lett. 168B (1986) 113. 

12. A.H. Mueller, Nucl. Phys. B348 (1991) 310. Nucl. Phys. B353 (1991) 44. 

13. T. Banks et al., Nucl. Phys. B347 (1990) 581. 

14. V.I. Zakharov, Nucl. Phys. B353 (1991) 683; 
G.Veneziano, Preprint CERN-TH.6399/92 (1992) 

15. M. Maggiore and M. Shifman, Nucl. Phys. B365 (1991) 161; Nucl. Phys. 
B371 (1992) 177. 

16. D. Diakonov and V. Petrov ,Non-Perturbative Isotropic Multiparticle Produc- 
tion in Yang-Mills Theory, Preprint RUB-TPII-52/93 (1993). 

17. V.A. Novikov, M.A. Shifman, A.I. Vainshtein and V.I. Zakharov Nucl. Phys. 
B174 (1980) 378. 

18. M.S. Dubovikov and A.V. Smilga, Nucl. Phys. B185 (1981) 109. 

19. I.I.Balitsky and V.M.Braun, Nucl. Phys. B311 (1988/89) 541. 

20. L.S. Brown, R.D. Carlitz, D.B. Creamer and C. Lee, Phys. Rev. D17 (1978) 
1583. 

21. ALEPH Collaboration, D. Buskulic et al., Phys. Lett. B307 (1993) 209. 



